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Chinese BeiDou Navigation Satellite System is officially operational as a regional constellation with five 
Geostationary Earth Orbit (GEO) satellites, five Inclined Geosynchronous Satellite Orbit (IGSO) satellites 
and four Medium Earth Orbit (MEO) satellites. Observations from the BeiDou Experimental Tracking 
Stations (BETS) and the IGS Multi-GNSS Experiment (MGEX) network from 1 January to 31 March 2013 
are processed for orbit determination of the BeiDou constellation. Various arc lengths and solar radiation 
pressure parameters are investigated. The reduced set of ECOM five-parameter model produces better 
performance than the full set of ECOM nine-parameter model for BeiDou IGSO and MEO. The orbit 
overlap for the middle days of 3-day arc solutions is better than 20 cm and 14 cm for IGSO and MEO in 
RMS, respectively. Satellite laser ranging residuals are better than 10 cm for both IGSO and MEO. For 
BeiDou GEO, the orbit overlap of several meters and satellite laser ranging residuals of several decimetres 
can be achieved. 



China is developing its own BeiDou Navigation Satellite System^ which is a useful constellation contributing 
to the GNSS community. BeiDou System development is divided into two phases, i.e. first the regional 
system and then the global system. Phase one has been accomplished with a constellation comprising five 
Geostationary Earth Orbit (GEO) satellites, five Inclined Geosynchronous Satellite Orbit (IGSO) satellites and 
four Medium Earth Orbit (MEO) satellites in orbit. After its experimental service in 2011, BeiDou has begun to 
provide official positioning, navigation and timing services for the Asia-Pacific area since December 2012. 
BeiDou satellites transmit navigation signals at three frequencies, i.e., 1561.098 MHz (Bl), 1207.140 MHz 
(B2) and 1268.520 MHz (B3). BeiDou aims at serving various fields including civilian and scientific applications, 
and thus high positioning quality is also required. Precise orbits are the foundation to expand its application and 
play an important role in complementing other GNSS. 

A number of studies have been carried out on precise orbit determination of BeiDou constellation. 
Montenbruck et al.^ presented initial assessment of the BeiDou regional navigation satellite system and the 3D 
orbit overlap accuracy was 1-10 m. In Ge et al.^, independent BeiDou orbit determination was carried out using 
three day arcs and nine solar radiation pressure (SRP) parameters of the ECOM model^. The 3D overlap of GEO 
and IGSO reaches 3.3 m and 0.5 m in RMS, respectively, and an almost constant disagreement in the along-track 
direction for GEO is found. Shi et al.^ adopted two-step GPS-assisted BeiDou orbit determination method, which 
is also used in orbit determination of GLONASS^, Galileo^ and QZSS^ and reached radial overlap precision at 
10 cm level for BeiDou GEO and IGSO satellites. Zhao et al.^ also conducted independent BeiDou orbit deter- 
mination using three day arcs, but estimated five SRP parameters of the ECOM model for all satellites and an 
additional constant acceleration in the along-track direction for GEO. The 3D overlap precision is about 1.8 m for 
GEO and 0.3 m for IGSO and MEO, and radial overlap precision is better than 0.1 m for all satellites. The SLR 
residuals, which mainly represent radial orbit error, reach about 0.7 m for GEO COl and 0.1 m for IGSO C08, 
which may indicate some unmodelled errors in orbit models, especially for GEO. Steigenberger et al.^° tested 
different arc lengths and SRP parameters in orbit determination of BeiDou GEO and IGSO satellites. In order to 
deal with strong correlations, only one SRP parameter in the direction of the Sun is estimated for the GEO 
satellites and three, five and nine parameters of ECOM SRP model are estimated respectively for the IGSO 
satellites. The results show that orbit consistency is on the several decimetre level for GEO and on the 1-2 dm 
level for IGSO. Independent SLR residual is about 0.15 m for GEO COl and 0.24 m for IGSO C08. The results also 
reveal that there might be problems with estimating only one direct SRP parameter for GEO because the orbit-fit 
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Figure 1 | Ground traces of BeiDou constellation and tracking stations. Ground traces of satellites are represented by red dots, blue and green curves for 
GEO, IGSO and MEO satellites, respectively. Stations from BETS are marked in brown and that from MGEX in blue. This figure is drawn 
using GMT software^^. 



RMS gets bigger with increasing elevations of the Sun above the orbit 
plane. Liu et al.^^ used one-step joint orbit determination of BeiDou 
and GPS as is recommended in Ineichen et al.^^ for GLOANSS and 
GPS. With ambiguity resolution for both GPS and BeiDou IGSO and 
MEO satellites, an overlap precision of about 1.5 m for BeiDou GEO 
satellites and 0.3 m for BeiDou IGSO and MEO satellites is reached. 
He et al.^^ conducted a series of experimental studies on the impact 
of the expanding tracking geometry, the involvement of MEO and 
integer ambiguity resolution in BeiDou orbit determination and 
indicated all these factors help to improve BeiDou orbit overlap 
precision. 

BeiDou Experimental Tracking Stations (BETS)^ is a continuous 
tracking network for BeiDou and GPS satellites established world- 
wide for scientific and engineering applications by the GNSS 
Research Centre of Wuhan University with cooperation partners 
since March 2011 and currently comprises fifteen stations. In Shi 
et al.^^, the analysis on geometry-based residuals of double difference 
phase and code shows that for BeiDou un-differenced carrier phase 
measurements the noise is about 2-4 mm on both Bl and B2, and 
less than 0.45 m for Bl CA code and about 0.35 m for B2 P code. The 
Multi-GNSS Experiment (MGEX)^' is organized by the IGS'" to 
track, collate and analyse all available GNSS signals including those 
from the Galileo, BeiDou and QZSS systems, as well as from mod- 
ernized GPS and GLONASS satellites and any space-based augmen- 
tation system (SBAS) of interest. Nine stations from the MGEX 
network are used in this study, which greatly extend the observation 
network for BeiDou satellites, especially for MEO satellites. The dis- 
tribution map of the selected stations from BETS and MGEX track- 
ing network for BeiDou orbit determination is shown in Figure 1. 



There is some difficulty and limitations in BeiDou orbit deter- 
mination. The tracking network is not well distributed worldwide 
like GPS and GLONASS. GEO satellites move slightly relative to 
ground tracking stations, which results in strong correlations and 
IGSO satellites move slowly compared to MEO satellites. The 
optimal precise force model of different types of satellites, especially 
SRP, still needs to be further investigated. In this study, we aim at 
precise orbit determination of BeiDou constellation based on BETS 
and MGEX network through various tests and assess the orbit accu- 
racy that can be achieved. The paper is organized as follows: data 
collection and processing method are described, then results about 
orbit overlap comparison and SLR validation are presented for dif- 
ferent tests, discussions are given in the last section. 

Data collection and processing. The observation data from the 
BETS and MGEX network from 1 January to 31 March 2013 are 
processed. Figure 1 shows the ground traces of BeiDou GEO, 
IGSO and MEO satellites and the tracking stations used in this 
paper. The station information is listed in Table 1. 

The Positioning and Navigation Data Analyst (PANDA) soft- 
ware has been developed at GNSS Research Centre, Wuhan 
University since 2000, as a multi-functional data processing software 
for various fields, such as precise orbit determination of GNSS and 
LEO satellites, gravity field recovery, real-time precise positioning 
services, and troposphere and ionosphere reconstruction. We have 
adapted the software package for BeiDou data processing in this 
study. There are generally three methods in precise orbit determina- 
tion for BeiDou Navigation Satellite System in multi-GNSS proces- 
sing. The first method is single-system mode, where all parameters 
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are estimated using only BeiDou observations. The second method is 
two-step GPS -assisted, where in the first step GPS observations are 
used to derive parameters common to both GNSS system (Uke sta- 
tion coordinates, receiver clock, troposphere delay parameters) 
which is then introduced as known quantities in BeiDou processing 
in the second step. The third method is one-step joint orbit deter- 
mination of BeiDou and other GNSS, like GPS, in one processing 
scheme. The second method can take advantage of the already avail- 
able precise GPS orbit and clock products^^ and result in a better 
estimation of BeiDou related parameters efficiently. Therefore, we 
adopt this method in this contribution. 

In the first step, GPS precise point positioning (PPP)^^ is carried 
out for the stations on the daily basis, using the ionosphere-free linear 
combination (LC) of LI and L2 code and phase observations with 
IGS final products, solving for daily station coordinates, epoch-wise 
receiver clock and 2-hourly zenith tropospheric delays (ZTD). These 
parameters are kept fixed in orbit determination in the second step 
using the ionosphere-free linear combination (LC) of Bl and B2 code 
and phase observations, estimating only BeiDou-related parameters: 
satellite initial position and velocity, SRP parameters, satellite clock 
corrections and differential code biases (DCB). These biases origin- 
ate from differences between the signals of BeiDou and GPS (B1/B2 
for BeiDou and L1/L2 for GPS) which cause inter-system biases 
(ISB) and inter-frequency biases (IFB). BeiDou and GPS may not 
have the same receiver antenna phase centre, but we can ignore its 
impact regarding current decimetre-level BeiDou orbit precision and 
millimetre-level phase centre diffrence^^. For all BeiDou satellites, 
since no accurate antenna phase centre (PCO) and phase centre 
variation (PCV) is available, we just use nominal PCO (0.634 m, 
— 0.003 m, L075 m) provided by the satellite manufacturer and 
ignore satellite PCV. Important options of processing strategy about 
observation model and force model are listed in Table 2. 

The differential code biases (DCB) have to be handled carefully in 
order to avoid biases in the final BeiDou orbits. In this paper, the 
DCB of the station Chengdu (CHDU) is fixed to zero as a reference 
and a piece-wise DCB with a resolution of one day is estimated for 
other stations. Mean DCBs and corresponding standard deviations 
for the experiment interval of three months are shown in Figure 2. 
The shaded stations are equipped with Unicore UR240 receivers 
while others are equipped with Trimble NetR9. The DCB consistency 



is about 15 ns for Unicore UR240 and 10 ns for Trimble NetR9. The 
mean standard deviation is about 1 ns for all stations which is stable 
in time for pseudorange noise level. 

The mean post-fit residuals of ionosphere-free linear combination 
of Bl and B2 code and phase observations from the BeiDou para- 
meter estimation are shown in Figure 3. The stations GRAC and 
BRST can only observe BeiDou IGSO and MEG satellites while the 
stations UNB3, ABMF and LMMF can only observe BeiDou MEG 
satellites. The mean phase residuals for all stations are about 1.8 cm, 
1.3 cm, and 1.4 cm for BeiDou GEO, IGSO and MEO satellites, 
respectively. The mean code residuals for all stations are about 
1.9 m, 1.8 m, and 1.9 m for BeiDou GEO, IGSO and MEO satellites, 
respectively. The GPS phase and code residuals of the first step are 
about 1 cm and 1 m respectively, i.e., the BeiDou residuals are larger 
than those in GPS, indicating some unmodelled errors in the BeiDou 



Table 2 I Observation models and force models 



Basic Observable 

Sampling rate 
Cutoff elevation 
Weighting 



Station phase center 

Satellite phase center 
Phase wind up 
Tropospheric delay 
Satellite clock error 
Receiver clock error 
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Tide displacement 
Relativity effect 
Earth gravity 
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N-body gravitation 
Atmospheric drag 
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7° 
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phase and code, respectively 
Elevation-dependent, 1 for E > 30°, otherwise 

2*sin(E) 
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same, igsOS.atx^*^ 
Only nominal PCO for BeiDou; igsOS.atx for GPS 
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GPT modeP\ GMF mapping function^^, 2-hourly 

White noise 

White noise 

One constant per day 

lERS Conventions 2003^^ FES2004^^ 

lERS Conventions 2003^^ 

EIGEN_GL04C^^upto 12x 12 
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Sun, moon and other planets 
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Figure 2 | DCB value and standard deviation. 



observations, like the satellite antenna phase centre and attitude 
behaviour. One has to note that it is generally difficult to directly 
reflect the orbit quality from observation residuals. 

Different arc lengths and sets of SRP parameters are tested for 
GEO, IGSO and MEO. The orbital arc lengths range from two to 
seven days, i.e., only one set of initial position and velocity and SRP 
parameters is estimated per arc. SRP is currently the largest 
non-conservative perturbation in orbit determination of navigation 



satellites which needs to be modelled accurately. The ECOM model^ 
comprises nine constants and periodic terms in three axes denoted as 
D, Y and B. D is along the satellite-sun direction and Y is along the 
spacecraft's solar panel axis. B is the direction perpendicular to D and 
Y and completes the right-handed coordinate system. SRP model has 
to take satellite attitudes into consideration. For BeiDou GEO satel- 
lites, orbit-normal attitude is always maintained^°. The spacecraft- 
fixed Z axis is pointing towards geo-centre and Y axis (along the solar 
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Figure 4 | Overlap comparison of orbit determination and prediction. 

panels) is oriented to the Sun and normal to the orbital plane. X axis 
is pointing to the along-track direction and completes the right- 
handed coordinate system. BeiDou IGSO and MEO satellites adopt 
yaw- steering attitude like GPS most of the time and only turn into 
orbit-normal attitude when the Sun elevation about the orbital plane 
is less than a threshold of about four degree^^'^^. The ECOM model is 
investigated for BeiDou satellites in both yaw-steering and orbit- 
normal attitude. In this contribution, ECOM 5 parameter and 
ECOM 9 parameter are tested for IGSO and MEO satellites while 
different subsets of ECOM model is tested for GEO with necessary 
modification. 

Results 

BeiDou IGSO and MEO. We assess BeiDou orbit quality in both 
internal consistency and external validation. For internal consis- 
tency, the comparison of satellite positions in the overlapped time 
span between different orbit solutions is used. For orbit arcs of 3, 5 
and 7 days, the last observed day and the first predicted day of orbit 
solutions with different arc lengths are compared with the middle 
day of the corresponding solutions. Figure 4 gives an example of 
overlap comparison for orbit solutions of 3 day arc. For orbit arcs 
of 2, 4 and 6 days, there are two days at the middle of each arc. The 



last observed day and the first predicted day of orbit solutions are 
compared with the first middle day of the corresponding solutions. 
For external validation, SLR observations are used to independently 
assess BeiDou orbit accuracy. The microwave-based orbits are kept 
fixed and residuals are computed between observed SLR values and 
the computed distance using BeiDou orbits. 

To evaluate the influence of the different number of ECOM para- 
meters on orbit internal consistency, the mean RMS values of overlap 
comparison in along-track, cross-track and radial direction for 
BeiDou IGSO are shown in Figure 5. We take the mean RMS values 
of the five BeiDou IGSO (C06-C10) satellites over the experiment 
time span of three months excluding orbit-normal attitude periods. 
The along-track direction shows the largest RMS values and the 
radial direction shows the smallest as GPS due to observation geo- 
metry. The results show that in general the RMS of ECOM 5 para- 
meter is smaller than that of ECOM 9 parameter. For ECOM 
5 -parameter model, the smallest RMS values are found at 3 day arc 
length and then show an almost linear increase with increasing arc 
length both in orbit determination and prediction. This indicates that 
3 day arc length is the best for ECOM 5 -parameter model in terms of 
internal consistency. For the ECOM 9 -parameter model, the smallest 
RMS values are found at 5 day arc length in radial direction for both 
orbit determination and prediction. The RMS values are found smal- 
lest at 3 day arc length for the cross -track direction and show an 
almost increasing trend in the along-track direction for orbit 
determination. 

As for BeiDou IGSO, the mean RMS values of overlap comparison 
in along-track, cross-track and radial direction for BeiDou MEO are 
shown in Figure 6. We take the mean RMS values of the four BeiDou 
MEO (C11-C14) satellites over the experiment time span of three 
months excluding orbit-normal attitude periods. The results also 
show that in general the RMS of ECOM 5 parameter is smaller than 
that of ECOM 9 parameter. For ECOM 5 parameter, the smallest 
RMS values are found at 3 day arc length in along-track and radial 
direction and then show an almost linear increase with increasing arc 
length in both orbit determination and prediction. For the cross-track 
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Figure 5 | (a) IGSO overlap comparison of the middle day with the last observed day. (b) IGSO overlap comparison of the middle day with the first 
predicted day. Orbit determination and one-day prediction with ECOM 5 parameter (red crosses) and ECOM 9 parameter (blue circles) are 
compared under different arc lengths. 
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Figure 6 | (a) MEO overlap comparison of the middle day with the last observed day. (b) MEO overlap comparison of the middle day with the first 
predicted day. Orbit determination and one-day prediction with ECOM 5 parameter (red crosses) and ECOM 9 parameter (blue circles) is 
compared under different arc lengths. 



direction, the smallest RMS values are found at 5 day arc length but 
are very close to that of 3 day arc length. We can draw the conclusion 
that 3 day arc length is the best for ECOM 5 parameter model for 
BeiDou MEO satelUtes in terms of internal consistency. For ECOM 9 
parameter, the RMS values show an almost linear decrease with 
increasing arc length both in orbit determination and prediction, 
but the values are close to each other when it comes to arc length 
of more than 5 days. The smallest RMS values in the radial direction 
are found at 5 day arc length in both orbit determination and 
prediction. 

From above analysis, for BeiDou IGSO and MEO satellites, the 
average RMS values show that ECOM 5 parameter model and 3 day 
arc length resulted in the best solution in terms of internal consist- 
ency in orbit determination and one- day prediction. In the next 
section, SLR data are analysed as an independent measurement to 
assess the BeiDou orbit. During the experiments, the BeiDou GEO 
COl, IGSO COS and CIO, MEO Cll equipped with laser retro- 
reflector arrays are being observed by several SLR stations coordi- 
nated by the International Laser Ranging Service (ILRS)^^. Outliers 



exceeding 0.3 m are excluded for IGSO/MEO and those exceeding 
2 m are excluded for GEO, and after removal of outliers there are 328 
normal points (NP) available for COl, 330 NPs for COS, 673 NPs for 
CIO and 430 NPs for Cll during the study period. The offsets and 
STDs of IGSO and MEO SLR residuals of the middle day for orbit 
solutions under different arc lengths and number of ECOM para- 
meters are given in Table 3. The result is generally in good agreement 
with the overlap comparison in radial direction in orbit determina- 
tion. For MEO Cll, the mean offsets and STDs are basically better for 
the ECOM 5 -parameter model than that for the ECOM 9 -parameter 
under different arc lengths. The smallest residuals are found at 3 day 
arc length for ECOM 5 parameter and at 5 day arc length for ECOM 9 
parameter, which is the same as that for the overlap comparison. For 
IGSO COS, the mean offsets and STDs are also better for ECOM 5 
parameter than that for ECOM 9 parameter at different arc length. 
The smallest residuals are found at 3 day arc length for ECOM 5 
parameter and at 5 day arc length for ECOM 9 parameter. For IGSO 
CIO, the overall RMS values are smaller for ECOM 5 parameter than 
that for ECOM 9 parameter under different arc lengths. The smallest 
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Figure 7 | (a) SLR residuals time series for COS. (b) SLR residuals time series for CIO. (c) SLR residuals time series for Cll. 



residuals are found at 3 day arc length for ECOM 5 parameter and at 
6 day arc length for ECOM 9 parameter. As an example, the times 
series of SLR residuals for the 3 -day orbits with the ECOM 5 -para- 
meter model are shown in Figure 7. The overall RMS values for COS, 
CIO, Cll are 8.1 cm, 8.9 cm, 7.3 cm, respectively, which are sUghtly 
larger than that of overlap comparison. 

From the analysis of orbit internal consistency (overlap compar- 
ison) and external validation (SLR residuals), the experiment shows 
that on average ECOM 5 parameter and 3 -day arc length result in 
better orbit for BeiDou IGSO and MEO satellites. Using this strategy, 
orbit overlap comparison of individual satellites in yaw- steering atti- 
tude mode in the experiment period is shown in Figure 8 (a). During 
yaw-steering period, the 3D precision is better than 20 cm and 14 cm 
for IGSO and MEO, respectively, and the radial precision is better 
than 5 cm and 3 cm for IGSO and MEO, respectively. During the 
experiments, several IGSO and MEO satellites went through orbit- 
normal attitude mode. The corresponding time periods are shown in 
Table 4. Using ECOM 5 parameter and 3 -day arc length strategy, 
orbit overlap comparison of these satellites in orbit-normal attitude 
mode is also obtained and shown in Figure 8 (b). The 3D precision is 
about 30 cm and 20 cm for IGSO and MEO, respectively, and the 
radial precision is better than 10 cm and 5 cm for IGSO and MEO, 
respectively. The results indicate that ECOM model with yaw- steer- 
ing can be applied to IGSO and MEO in orbit-normal attitude mode 
although the orbit precision is slightly worse than that in yaw- steer- 
ing mode, which may indicate deficiencies related to the attitude 
modelling. 



BeiDou GEO. Due to minimal changes of the observation geometry, 
strong correlations occur among the SRP parameters, orbital 
elements, satellite clocks and ambiguities in BeiDou GEO orbit 
determination. To deal with these correlations in GEO orbit 
determination, Steigenberger et al.^° estimated only one direct 
ECOM parameter in the direction of the Sun. The internal 
consistency is on the several decimetre level for GEO satellites and 
the overall RMS of SLR residuals is 14.7 cm from 10 normal points. 
The results also reveal that estimating only one direct SRP parameter 
might be a too simple model because the orbit-fit RMS gets bigger 
with larger elevation of the Sun above the orbit plane. Orbit-normal 
attitude is always adopted for BeiDou GEO and the solar panel is not 
perpendicular to the incidence sunlight. Zhao et al.^ estimated 
ECOM 5 parameter and an empirical constant acceleration in the 
along-track direction to compensate for the insufficiency of the 
ECOM for GEO. The overall overlap comparison is better than 
200 cm in 3D and better than 10 cm in the radial direction and 
the overall RMS of SLR residuals is 68.5 cm from 25 normal 
points. The study also shows that the significant mean deviation 
from zero in the SLR residuals may be problematic. We propose 
for GEO a modified ECOM 5 parameter model with a virtual Y 
axis, resulting in three orthogonal axes denoted as D, Yv and B. In 
ECOM model for BeiDou GEO, the axis Y is perpendicular to the axis 
B, but it is not perpendicular to the axis D due to the orbit-normal 
attitude if the sun elevation about the satellite orbital plane is not 
zero. In our modified ECOM model, the direction of the virtual axis 
Yv is the cross product of the ECOM axes B and D. With this 
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Figure 8 | (a) overlap comparison in yaw- steering mode, (b) overlap comparison in orbit-normal mode. 
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Table 4 


Orbit-normal attitude period during the study 


PRN 


DOYof2013 


C07 


28-41 


COS 


1-9 


CIO 


27-40 


Cll 


16-24 


C12 


16-24 



modification, the virtual axis Yv has an angle with respect to the solar 
panel axis which changes slowly with the elevation of the Sun above 
the GEO orbit plane. The virtual axis Yv is thus able to absorb the 
SRP force perpendicular to the solar panel introduced by reflected 
and diffusely scattered light photons. These three SRP models used in 
Steigenberger et al.^°, Zhao et al.^ and this paper are denoted as model 
1, model 2 and model 3, respectively. With orbit maneuvers^° 
excluded (three for COl, C02, C03, COS respectively and one for 
C04), these three SRP models are compared in GEO orbit 
determination with 3 -day arc length and average overlap 
comparison for all GEOs and SLR residuals for COl are given in 
Table 5 and Figure 9, respectively, and the results of other arc 
length show similar patterns. 

There are 328 normal points (NP) available for COl during the 
experiment period. The COl satellite is observed by Changchun sta- 
tion, Shanghai station in China and Yarragadee station in Australia. 
The relation between expected SLR residual and orbital error in 
along-track, cross-track, radial direction can be derived. The orbital 

error vector e is expressed as follows, where A, C, R are orbital error 

values and a , c , r areunitvectorsinalong-track, cross-track, radial 
direction respectively. 

1 =A'a ^Cl ^Rr (1) 

The unit vector in the line-of-sight direction from the GEO to one 

station is denoted as / . The projection of the orbital error vector 
onto the line-of-sight direction can be obtained by the dot product as 
follows. 

e .7 =(A a + C c r ).7 = A a-T + C c - T r - T (2) 

Since the GEO almost does not move relative to ground stations, a , 

c , r and / change little and can be computed using approximate 
coordinates. We derive the approximate formula for the three sta- 
tions as follows, which shows that SLR residuals mainly reflect the 
orbital error in radial direction. 

CHAL = 0.03A - 0.14C - 0.99i^ 

SHAl=0.05A-0.11C-0.99i^ (3) 

YAKL = 0.06A + 0.06C - 0.99i^ 

As far as internal consistency is concerned. Model 1 has the largest 
overlap difference in radial direction which is about one order of 
magnitude larger than that of Model 2 and Model 3. The 3D pre- 
cision is about 4 m, 1.5 m and 1.8 m on average, respectively. The 



Table 5 


SRP model comparison 


for BeiDou GEO 






Overlap 


comparison 


(cm) 


SLR residuals 


(cm) 


Along 


Cross Radial 


Mean STD 


RMS 


Model 1 


387.6 


21.2 


88.9 


-31.9 40.8 


51.8 


Model 2 


146.4 


5.8 


6.2 


-70.2 33.9 


77.9 


Model 3 


173.7 


46.3 


9.2 


-57.2 32.6 


65.8 



biggest error occurs in along-track direction and it is interesting to 
note that precision is better in cross-track direction than that in 
radial direction for Model 1 and Model 2. When it comes to SLR 
residuals. Model 1 shows the smallest mean and RMS value while 
Model 2 shows the largest mean and RMS value, which is contra- 
dictory to the conclusion from the overlap comparison in radial 
direction. For Model 3, STD of SLR residuals is the smallest and 
RMS value is 14 cm larger than that of Model 1. The discrepancy 
between overlap comparison and SLR residuals exists in these three 
models, and the least pronounced is for Model 1. This may indicate 
some unmodelled errors related to orbit dynamics or spacecraft. 
There should be a balance in the number and composition of SRP 
parameters for BeiDou GEO to account for both strong correlation 
and accurate description of orbit dynamics. This is still an open 
research problem and needs further investigations. Other SRP model 
and possible variations should also be studied in BeiDou orbit deter- 
mination, like adjustable box-wing modeP^ and GSPM modeP^. 
Stable satellite clock can be used as indicators for orbit modelling 
and satellite clock model constraint can improve orbit accuracy^^. 
Clock modelling may improve BeiDou orbit accuracy with proper 
constraints, especially for BeiDou GEO satellites. Currently the 
BeiDou satellite clocks are not stable enough^'^° and the future clocks 
are expected to do better. The proper constraints for the BeiDou 
satellite clock model should be adjusted accordingly and weak con- 
straints should probably be applied to the BeiDou satellites which 
need further investigations. Due to the fact that GEO is basically 
stationary relative to the regional observing stations, the satellite 
orbit and clock is highly correlated. Carolipio et al.^^ and Grewal et 
al.^^ proposed an effective technique to decouple the clock error as 
known quantity and explored the great improvement in GEO orbit 
determination. As for BeiDou satellite, independent operational time 
synchronization system is realized by comparing two-way up -link 
and down-link pseudo-range. It is expected that these observations 
are of high quality, and if the data are available, it can also be used to 
derive satellite clock as known quantity to improve BeiDou orbit 
determination, especially for GEO. 

Discussion 

We processed three months of observations from 23 stations in BETS 
and MGEX network for precise orbit determination of BeiDou 
regional navigation constellation of 5 GEO, 5 IGSO and 4 MEO 
satellites. Two-step GPS-assisted BeiDou orbit determination 
approach is implemented in the PANDA software and different 
arc lengths and sets of SRP parameters are tested for GEO, IGSO 
and MEO. Internal consistency and SLR validations show that for 
IGSO and MEO the orbit can reach accuracy on the level of 1-2 
decimetres and the reduced set of ECOM 5 parameters model 
resulted in a better performance than the full set of ECOM 9 para- 
meters model. For BeiDou GEO satellites, the tests show that internal 
precision of several meters and satellite laser ranging residuals of 
several decimetres can be achieved. There are still some problems 
in BeiDou GEO orbit determination and the discrepancy between 
overlap comparison and SLR residuals indicate some unknown sys- 
tematic errors, such as satellite attitude maintenance deviation. 
Further efforts may involve two aspects. One is reducing the strong 
correlation among parameters due to poor observation geometry and 
possible solutions are satellite clock modelling or determination via 
extra observation, like independent time synchronization system. 
The other is to try to model the orbit dynamics accurately and other 
empirical or analytical SRP models may be adapted for BeiDou orbit 
determination. The proper and better SRP model for BeiDou IGSO 
and MEO satellites should also be further investigated. 

As BeiDou is still under development, a larger tracking network 
with more stations and more even distribution around the world can 
provide better coverage for all satellites, especially for MEO. The global 
network can help to improve BeiDou service, like the estimation and 
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Figure 9 | (a) SLR time series for GEO COl with SRP model 1. (b) SLR time series for GEO COl with SRP model 2. (c) SLR time series for GEO COl with 
SRP model 3. 



refinement of the PCO and PCV model of satellites. Moreover, the 
three-frequency observation data should be fully exploited. Similarly, if 
more BeiDou satellites can be tracked by more SLR stations, it can 
provide more valuable data and powerful tool to analyse and improve 
BeiDou precise orbit determination. 
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